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Abstract 

We consider the dynamics of the scalar and pseudoscalar glueballs in the Quark-Gluon 
Plasma (QGP). By using the instanton model for the QCD vacuum we give the arguments 
that the nonperturbative gluon-gluon interaction is qualitatively different in the confine¬ 
ment and deconfinement phases. Based on this observation it is shown that above T c the 
values of the scalar and pseudoscalar glueball masses might be very small. The estimation 
of the temperature of scale invariance restoration, at which the scalar glueball becomes 
massless, is given. We also discuss the Bose-Einstein condensation of the glueballs and 
the superfluidity of the glueball matter in QGP. 
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1 Introduction 


In spite of the tremendous experimental and theoretical efforts in the investigation of the 
properties of quark-gluon plasma (QGP) at large temperatures and densities, it is not 
clear so far what is the fundamental QCD mechanism leading to the unusual behaviour 
of the matter produced in heavy ion collisions at high energies at the RHIC and LHC. 
Indeed, there is strong evidence that even above T c non-perturbative QCD effects are very 
important and, in particular, they are responsible for the phenomenon of the so-called 
strongly-interacting Quark-Gluon Plasma (sQGP) [Tj observed at the RHIC and LHC. 
One of the important outputs of the Lattice QCD calculation at finite T is the fact that the 
values of the deconfinement and chiral restoration temperatures are approximately equal. 
That means that above T c one cannot expect the pion to be the Goldstone boson with 
zero mass in the chiral limit. Instead, a rather massive (M n cs 2 M q (T c ) ~ 6 T c ~ 1 GeV) 
pion state appears above T c m 150 MeV. Therefore, this lowest mass quark-antiquark 
state can not give a significant contribution to the Equation of State (EoS) of the QGP. 

One of the fundamental issues of QCD as the theory of strong interactions is the 
understanding of the role of gluonic degrees of freedom in the confinement and deconfine¬ 
ment regimes. Thus, it is known that at low temperatures gluons play an important role 
not only in the dynamics of usual hadrons. In particular, they can form bound states 
called glueballs (see review 0). However, one cannot expect a significant contribution of 
glueballs at T <T C to EoS of the hadron gas due to their large masses M G 3> T. On the 
other hand, it was found by lattice calculation that unlike the quark condensate, the 
gluon condensate does not vanish at T > T c . We should stress that this condensate plays 
a fundamental role in the formation of the bound glueball states. This role is similar to 
the role of the quark condensate in the appearance of the massive mesons and baryons 
built of the light u, d, and s quarks. Therefore, a non-zero value of the gluon condensate 
above T c is a strong signal of the existence of the glueballs in the deconfinement phase. 

However, it is evident that the properties of the glueballs, in particular, their masses 
and sizes should change in the QGP due to the temperature dependence of the gluon 
condensate OH and the change of the topological structure of the QCD vacuum at 
T > T c \ 5], 1]. The attempts to estimate the value of the glueball masses at finite T were 
done in the different models in the papers PEI US CZl ESI IMj and in the lattice |2T, [22] , 

In this Letter, we consider the lowest scalar and pseudoscalar glueball states at finite 
T in the QGP environment within the effective model based on the instanton picture for 
the QCD vacuum. It is shown that their masses strongly decrease above T c and become 
very small at the temperature of the scale invariance restoration T sca i e 1 GeV. The 
possibility of the Bose-Einstein condensation and superfluidity of the glueball matter is 
under discussion. 
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2 Nonperturbative quark-quark, quark-gluon and gluon- 
gluon interactions below and above T c induced by 
instantons 


It is expected that the instantons, a strong vacuum fluctuation of gluon fields, play a very 
important role in the dynamics of the glueballs below T c (see review [23]). The effective 
interaction induced by the instanton between quarks at T = 0 is well known. This is a 
famous t’Hooft interaction which for Nf = 3 (Fig. la) and N c = 3 is given by the formula 
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where, mf ir is the quark current mass, q R) L = (1 ± r y 5 )q(x)/2, = qi R X a qiL, ji^ v = 

<jiR&[iu^ a qiL, P is the instanton size and n(p) is the density of instantons. For Nf — 2 the 
t’Hooft interaction for zero current quark mass is much simpler 
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This Lagrangian can be obtained from Eq. (JT]) by connecting the strange quark legs through 
the quark condensate. The two Lagrangians Eq. (]Tj) and Eq.Q are considered as the bases 
for different versions of the Nambu-Jona-Lasinio (NJL) models [24] , 

The effective quark-gluon interaction induced by the instantons is (see, for example 

mm) 
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1 In this section, all Lagrangians are written in the Euclidian space-time. 


(3) 
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where U is the orientation matrix of the instanton in the SU( 3) c color space. From this 
Lagrangian one can obtain the interaction between quarks and gluons which contributes 
to the mixing between scalar and pseudoscalar glueballs, quarkoniums and tetraquark 
states 
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where 
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where G^ v (x) = \/2e llva pG a a ^[x\ and Cfj and L^a are the t’Hooft interaction induced 
by the instanton and antiinstanton, respectively. 



Figure 1: The quark-quark a), quark-gluon b), c) and gluon-gluon interaction d) induced 
by the instanton for the Nf = 2 case. The symbol x means connection through the quark 
condensate or by the currect quark mass. 


The interaction between gluons in the scalar and pseudoscalar channels induced by 
the instanton is 


7T 6 /• ; n(p)p 8 

99 ~ 80 J P «2(p) 
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Due to factor 2 r n 2 / g s in Eq.(j3j) each two gluon legs produce a large enhancement factor 
7r 6 /«^ in the gluon-gluon interaction induced by instantons. So one can expect that this 
type of the interaction is much stronger in comparison with the quark-quark case. 

We would like to emphasize that the instanton induced interaction is very sensitive to 
the parity of the glueball and quarkonium states and, in particular, it leads to the mass 
splitting between scalar and pseudoscalar glueballs [3D] 07] • Indeed, the single instanton 
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contribution to the difference of two correlators of glueball currents with opposite parities 
is given by 


AfI(Q 2 ) G = i J d 4 xe iqx (< 0|T0 5 (o;)Os(0)|0 > - < 0\T\0 P (x)0 P (0)\0 >) 

= 2 6 vr 2 J dpn(p)(pQ) 4 K 2 2 (pQ), (7) 

where the glueball currents for the scalar and pseudoscalar states are the following: 

Os(x) = asG^WG^x), (8) 

0 P (x ) = a s G^(x)G; v (x). (9) 



Figure 2: The quark-quark a) and gluon-gluon interaction b) induced by instanton- 
antiinstanton molecule. 

At T < T c the main contribution to the quark-quark and gluon-gluon interaction 
comes from disordered instantons, Fig.l, in the so-called ’’random” phase of the instanton 
liquid. The phase provides spontaneous chiral symmetry breaking in QCD (see reviews 
[231124] ). In this case, the chiral symmetric contribution, which is related to the correlated 
instanton-antiinstanton molecules, Fig.2, is expected to be small due to the small packing 
fraction of instantons in the QCD vacuum [23]. Above T c , where the chiral and 1 / 4 ( 1 ) 
symmetries are restored, the situation is opposite. Indeed, in this case the contribution 
from the chirality violated interaction, Fig.l, is proportional to the product of the current 
quark masses and should be small. I 11 the chiral symmetric phase, the leading contribution 
comes from the strongly correlated instanton-antiinstanton molecules, Fig.2, as supported 
by the calculations in [2B]. The contribution from the eight-quark interaction presented 
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b) 


Figure 3: The eight-quark a) and eight-gluon interactions b) induced by instanton- 
antiinstanton molecule for the Nf — 3 case. 

in Fig.3a does not give a significant contribution to the hadron spectroscopy but it is 
important for the stability of the vacuum of the NJL model with the instanton induced 
interaction mm- 

The effects of quark-quark interaction induced by instanton-antiinstanton molecules 
near T c were studied in [32] 133]. ft was shown that they might give a significant contri¬ 
bution to the binding energy of the quark-antiquark states in the deconfinement phase. 
The main effect, in comparison with the zero temperature case, comes from the strong 
polarization of instanton-antiinstanton molecules in the time direction at the high tem¬ 
perature. 

3 Lowest mass scalar glueball above T c 

At T = 0 there are several scalar meson states /o(600), /o(980), /o(1500) and /o(1710) 
with a possible admixture of gluons pj^. Unfortunately, even the lattice calculation, 
which is based on the first principles of QCD, cannot give an accurate prediction for 
the scalar glueball masses with unquenched quarks. The main problem here is in the 
large contribution coming from the so-called disconnected diagrams [10] . In particular, 
they might also be responsible for the strong quarkonium-glueball mixing. Within the 
instanton model such mixing was shown to be large in the different approaches in ram 
Therefore, it is quite difficult, even impossible, to find a pure glueball state at T = 0. 
Several arguments based on the analysis of sigma-meson contribution to different reactions 

2 We do not include in this list the /o(1370) state due to its controversial situation in the experiment 
(see the recent discussion in j5] ). 
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were given to consider the lowest mass scalar sigma-meson / 0 (600) state as the state with 
a large mixture of the glueball umm- In the line of these studies, we will also assume 
that at T = 0 the lightest glueball state is the sigma meson state with the mass m a ~ 450 
MeV. It is also well known that the masses of the hadron states are changing only a 
little in the interval 0 < T < T c due to a small change of the values of quark and gluon 
condensates in this region. Therefore, the mass of sigma should be around 450 MeV at 
T f« T c . It has been discussed above that in the chiral symmetric phase above T c the main 
contribution to the gluon-gluon interaction is related to the formation of the instanton- 
antiinstanton molecules. The example of these interactions is presented in Fig.2b and 
Fig.3b. Similarly to the case of the NJL model with the eight-quark interaction based 
on the instanton induced multiquark interaction [301 3T], the effective Lagrangian in the 
gluon sector has the form 

£(*) = - V(4>), (10) 

where 

r„( 4 ) =(11) 

and $ is now the scalar glueball held. 

Due to the scale anomaly in QCD, the trace energy momentum tensor is not zero but 
it is proportional to the value of the gluon condensate 

= --< —G 2 > . ( 12 ) 

M 8 7r v ; 

By using the parton-hadron duality principle we obtain the relation between the gluon 
condensate in QCD and the minimum of the potential, EqlTTl in the glueball sector 

T% = 4V(<F 0 ). 

To estimate the coupling A, we use the relation 

27rr?7.Q 

9 < a s G 2 > ’ 

and the value of the gluon condensate m 

< a s G 2 0.076V1' 4 . 

For the mass of the glueball at T « T c m 0 ~ 450 MeV we have A = 0.41. Now we are in 
a position to calculate the temperature dependence of the scalar glueball mass above T c . 
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At the finite temperature, the effective potential in the one loop approximation for the 
A$ 4 theory is given (see for example, US!) 
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Figure 4: The ratio of m\(T))/ml as a function of T/T c . 
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For the trace anomaly we obtain 


m%{T) 

4A 


= 4 = -A $i in (T) = 


The temperature dependence of the gluon condensate in this model is given by the formula 


< — G 2 (T) >= 
7r 


2mj(T) 

9A 


(18) 


In FigHl the temperature dependence of the mass of the scalar glueball is presented. So 
one can see that at T sca i e « 6T C 0.9 GeV the mass of the scalar glueball vanishes and 
the scale invariance is restored. This can be treated as the appearance of the massless 
dilaton held in the limit —> 0 (see the recent discussion of the dilaton in [43]). 

We should emphasize that the justification of the use of the point-like effective inter¬ 
action in Eq. flTTT) can be related to a very small size of the scalar glueball in the QGP. 
Even at T = 0, due to the strong attraction between gluons induced by the instantons, 
the size of the scalar glueball is very small, R g & 2/3 p c ~ 0.2 fm, where p c ~ 0.3 fin is the 
average instanton size in the QCD vacuum [33]. This small size was also confirmed in the 
lattice calculation [35], [36]• Above T c a similar phenomenon happens as well. However, 
in this case the instanton-antiinstanton molecules produce a very strong attraction in the 
scalar glueball channel. Therefore, we expect that the size of the scalar glueball should 
be smaller than the size of instantons in the QGP, which is cutting at 


Rg « P\T) « 


1 

3tt 2 T 2 ’ 


(19) 


at finite T and N c = Nf = 3 [39]. We would like to emphasize that this size is smaller 
than the perturbative Debye screening length in the QGP 


Apr) = 1/MpT) > -D (20) 

where a s (p c ) ~ 0.5 [23] for T = 0 and M|,(T) = g 2 s (N c /3 + Nf/6)T 2 [38] was used. 
Therefore, we come to the important conclusion that the Debye screening cannot destroy 
the binding of the scalar glueball in QGP. 


4 Pseudoscalar glueball above T c 

There are also several candidates for the pseudoscalar glueball below 2 GeV at T = 0: 
?7(1405),7/(1475),7/(1760) and A"(1835). In the lattice quenched calculation the lowest 
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mass pseudoscalar glueball was predicted to have the mass M 0 -+ ~ 2.6 GeV Bound The 
large difference between scalar and pseudoscalar glueball masses m 0 -+ (T = 0) —m 0 ++ (T = 
0) ~ 1 GeV at low temperature can be explained by the large single instanton contribution 
to these channels (see discussion in [23]). Indeed, it gives a strong attraction in the scalar 
channel and a strong repulsion in the pseudoscalar state. After restoration of the chiral 
symmetry the masses of scalar and pseudoscalar glueballs should be equal in the limit of 
the zero light quark masses, m u = m d = m s = 0. In this case, instanton-antiinstanton 
molecules give the same strength of attraction in the both states. Therefore, the leading 
contribution to the splitting between the masses of two states at T > T c is determined by 
the density of the single instantons, which is proportional to the product of the current 
masses of the light quarks. We should mention that even at large temperature the single 
instanton contribution is repulsive in the pseudoscalar glueball case, and its collapse to the 
massless state is not allowed. This is a very important difference between the properties 
of pseudoscalar and scalar glueballs at the temperature of the scale invariance restoration. 

We can estimate the mass of the pseudoscalar glueball at T = T sca i e by using the 
instanton model. The difference between the scalar and pseudoscalar glueball masses 
above T c is determined by the instanton density n(p(T )) 

n{p{T)) ~ m u m d rn s (p(T)j , (21) 


where b 0 = HN c /3 ■ 
have 

^ 0 - + ( T SC ale ) 


2Nf/3 and rn l are the current masses of the quarks. Finally, we 


m u m d m s / 1 

m* u m* d m* V ^3 TrT scale p, 


(m 0 -+(T = 0) - m 0 ++ (T = 0)), (22) 


where m* = , rn l + m* is the so-called effective mass of the quark in the instanton vacuum 
related to the quark condensate. The use of m u = m d ~ 4.5 MeV, m s ~ 100 MeV [42], 
m* = 170 MeV [23], p c = 1/600 MeV -1 and m 0 -+(T = 0) — m 0 ++(T — 0) « 1 GeV gives 
the estimation 


m 0 -+(T sca ie) ~ 0.1 eV. (23) 

Therefore, the mass of the pseudoscalar glueball at T = T sca i e is finite but very small. 

3 The difference between the observed and lattice masses might be related to the strong mixing between 
the glueball and quarkonium states [T2] ■ 
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5 Bose-Einstein condensation of the light glueballs 
and superfluidity of glueball matter in QGP 

The Bose-Einstein condensation (BEC) of the identical bosons is the well-known phe¬ 
nomenon in both theoretical and experimental physics. Recently, it was suggested that 
the over-occupied initial state of gluons created during relativistic heavy ion collisions 
might lead to the BEC of gluons HU ES SSI HZ]- However, gluons have the color charge 
and spin and, additionaly, the number of gluons is not conserved. All of these features, 
lead to a very complicated study of the possible formation of the gluon BEC in the QGP. 
On the other hand, the light glueballs with zero spin and color can easily form BEC at 
rather high temperature because for the given boson number density n the BEC temper¬ 
ature depends on the mass as [48] 


n 2/3 

T bec « 3.31-, (24) 

m 

for the non-relativistic case. For the arbitrary boson mass the critical BEC density is 
related to the temperature by the equation (see, for example, [49] and references therein) 


riBEC = 


d 3 k 1 

(27t) 3 exp((\Zk 2 + m 2 — m)/T B Ec) ~ 1 


(25) 


In FigE we show the critical glueball density of the scalar and pseudoscalar glueballs 
in the QGP by using the T-dependent glueball mass from Eq. dlGl) 0. This density can 
be compared with the density of QGP at the time r produced at the RHIC, which was 
estimated by using the multiplicity of the meson production in the Au-Au central collision 

PI 

13 

riRHic « — ( fm ~ 3 ), (26) 

r 

where r is in fm. For LHC Pb-Pb central collisions at y/s^N = 2.76 GeV the multiplicity 
is approximately twice larger. Therefore, the density is 


riLHC ~ — (fm 3 ). 

T 


(27) 


For QGP in the thermal equilibrium we can estimate the contribution of the glueballs to 
the total density. In this case, we have approximately 16 gluon and 24 quark degrees of 
freedom (Nf = 2) and 2 for glueballs. For the gluons and quarks with the mass M qg & 3 T 

4 We neglect a tiny difference between the values of the scalar and pseudoscalar glueballs in the QGP. 
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n(T) 



Figure 5: The T-dependency of the BEC density of glueballs in the 1 /fm 3 units. 

[51] . we have for scalar (pseudoscalar) glueball contribution to the total density at the 

RH1C 

4»w « ;(/™- 3 ), (28) 

for the initial temperature T^ HIC f« 300 MeV. 

For the LHC energy ( Tf HC fa 400 MeV) the estimation is 

« ;(/>"- 3 )- (29) 

By using the Bjorken model for the QGP expansion [52] 

tT 3 = r 0 T 0 3 , (30) 

with the thermalization time To fa 1 fm, we can estimate BEC temperature for the RH1C 
and LHC 

Trec C « 200MeV, Tj(a 270 MeV. (31) 

Therefore, the formation of the glueball BEC is possible at both RHIC and LHC heavy ion 
collisions. We should also mention that there is a large gap between the mass of the light 
glueball and the mass of the excited state in the system of two gluons in the QGP, which 
is m gg fa 6 T > 1 GeV. It is well know that such a gap should lead to the phenomenon 
of the superfluidity of the BEC matter. So we arrive at the conclusion that the QGP at 
the RHIC and LHC might be considered as the mixture of three matters. One of them is 
the usual ’’normal” QGP matter consisting of quarks and gluons and other ones are two 
superfluid matters of very light scalar and pseudoscalar glueballs. 
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Figure 6: The production of the scalar S and pseudoscalar P glueballs in the nonequi¬ 
librium QGP. 


There is also an additional mechanism which can lead to the abundance of the light 
glueballs production in ultra-relativistic heavy ion collisions. Indeed, at very high temper¬ 
atures at the initial stage of the production of the quark-gluon matter, the pair production 
of glueballs is possible by the fusion of two gluons, Fig.6. The effective interaction re¬ 
sponsible for such production is 

£ = A 2 (T)a a G%,G%,(S 2 + P 2 ), (32) 

where S(P) are the scalar (pseudoscalar) fields. In the mean field approximation we have 

A JT) ~ ms ' p(r)2 ~ 1 

21 ' c,,G; u G“„(T) m|_ p (T) 

Therefore, one might expect a strong enhancement of the ultra-light glueball production 
in the phase where the quark-gluon plasma is far away from the equilibrium. 


(33) 


6 Conclusion 

In summary, we considered the properties of scalar and pseudoscalar glueballs in the 
Quark-Gluon Plasma created in relativistic heavy-ion collisions. Based on the instanton 
model for the QCD vacuum we gave the arguments in favor of the existence of very 
light scalar and pseudoscalar glueball states above the temperature of the deconfinement 
transition. The estimation of the temperature of the scale invariance restoration, at which 
the scalar glueball becomes massless, was given. We also discussed the mechanism of the 
Bose-Einstein condensation and the superfluidity of the scalar and pseudoscalar glueball 
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matter in the QGP. We also showed the possibility of the abundance glueball production 
at the initial stage of the QGP formation. The influence of this phenomenon on the fast 
thermalization of the QGP is the subject of our future investigation. 
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